cpi uiaii 



Disorder correlations at smeared phase transitions 

Christopher Svoboda, David Nozadze, Fawaz Hrahsheh and Thomas Vojta 
Department of Physics, Missouri University of Science and Technology, Rolla, MO 65409, USA 



PACS 05 . 30 . Rt - Quantum phase transitions 
PACS 75 . 10 . Lp - Band and itinerant models 
PACS 75. 10. Nr - Spin-glass and other random models 

Abstract - We investigate the influence of spatial disorder correlations on smeared phase tran- 
sitions, taking the magnetic quantum phase transition in an itinerant magnet as an example. We 
find that even short-range correlations can have a dramatic effect and qualitatively change the 
behavior of observable quantities compared to the uncorrelated case. This is in marked contrast 
to conventional critical points, at which short-range correlated disorder and uncorrelated disorder 
lead to the same critical behavior. We develop an optimal fluctuation theory of the quantum phase 
transition in the presence of correlated disorder, and we illustrate the results by computer simu- 
lations. As an experimental application, we discuss the ferromagnetic quantum phase transition 
in Sri-^Ca^RuOs. 



Introduction. — Quenched disorder has various im- 
portant consequences in condensed matter. For example, 
disorder can change the universality class of a critical point 

• [HH] or even change the order of a phase transition [3H5] . 
In theoretical studies, the disorder is often assumed to 

. be uncorrelated in space even though many sample prepa- 
ration techniques will produce some degree of correlations 
between the impurities and defects. As long as the cor- 
relations are short-ranged, i.e., characterized by a finite 

. correlation length ^is, this assumption is usually justified 

' if one is interested in the universal properties of critical 
points. (There are exceptions for special, fine-tuned lo- 
cal correlations [6]). The reason why short-range corre- 
lated disorder leads to the same behavior as uncorrelated 
disorder can be easily understood within the rcnormaliza- 
tion group framework. Under repeated coarse graining, a 
nonzero disorder correlation length £dis decreases without 
limit. The disorder thus becomes effectively uncorrelated 
on the large length scales that determine the critical be- 
havior. 

A formal version of this argument follows from the Har- 
ris criterion [7]. It states that a clean critical point is 
stable against weak uncorrelated disorder if its correlation 
length critical exponent v fulfills the inequality dv > 2 
where d is the space dimensionality. If the inequality is 
violated, the disorder is relevant and changes the critical 
behavior. According to Weinrib and Halperin [5], spa- 
tially correlated disorder leads to the same inequality as 
long as its correlations decay faster than r~ d with distance 



r. Thus, short-range correlated disorder and uncorrelated 
disorder have the same effect on the stability of a clean 
critical point. 

In this Letter, we demonstrate that spatial disorder 
correlations are much more important at smeared phase 
transitions, a broad class of classical and quantum phase 
transitions characterized by a gradual, spatially inhomogc- 
ncous onset of the ordered phase [9] . Specifically, we show 
that short-range correlated disorder and uncorrelated dis- 
order lead to qualitatively different behaviors. The dis- 
order correlations do not only influence quantities usually 
considered non-universal such as the location of the phase 
boundary, they also change the functional dependence of 
the order parameter and other quantities on the tuning 
parameters of the transition, as indicated in Fig. [TJ We 
propose that this mechanism may be responsible for the 
unusually wide variations reported in the literature on the 
properties of the ferromagnetic quantum phase transition 
(QPT) in Sri-a.Ca.RuOa. 

In the following, we sketch the derivation of our theory, 
compute obscrvables, and illustrate them by simulations. 
We also discuss the generality of our findings, and we com- 
pare them to experiment. 

Smeared quantum phase transition. For defi- 
niteness, we consider a magnetic QPT in a metallic sys- 
tem with Ising order parameter symmetry. In the absence 
of quenched disorder, the Landau-Ginzburg- Wilson free 
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Fig. 1: (Color online) Schematic of the zero-temperature 
magnetization-composition curve (M vs a;) and the finite- 
temperature phase boundary (T c vs a;) at a smeared quan- 
tum phase transition in a random binary alloy A\- X B X . The 
cases of uncorrelated, correlated, and anti-correlated disorder 
are contrasted. 

energy functional of this transition is given by [TUlfTT) 

S = J dydz ^(y)T(y, z)^(z) + u J dy i'\y) , (1) 

where ip is the order parameter field, y = (y, r) comprises 
(i-dimcnsional spatial position y and imaginary time r, 
the integration means J dy = J dy Jq dr. and u is the 
standard quartic coefficient. The Fourier transform of the 
Gaussian vertex T(y, z) reads 

r(q,w n ) = r + £oq 2 +7o(q)K| . (2) 

Here, r is the distance from criticality0 £o is & microscopic 
length, and uj n is a Matsubara frequency. The dynamical 
part of r(q, uj n ) is proportional to \ui n \. This reflects the 
Landau damping of the order parameter fluctuations by 
gapless electronic excitations in a metallic system. The 
coefficient 70 (q) is q-independent for an antiferromagnetic 
transition but proportional to l/|q| or l/|q| 2 for ballistic 
and diffusive ferromagnets, respectively. 

We now consider a random binary alloy Ai_ x B x con- 
sisting of two materials A and B. Pure substance B has 
a non-magnetic ground-state, implying a positive distance 
from quantum criticality, re > 0. Substance A has a mag- 
netically ordered ground state with r\ < 0. By randomly 
substituting B atoms for A atoms, one can drive the sys- 
tem through a QPT from a magnetic to a nonmagnetic 
ground state. 

Due to statistical fluctuations, the distribution of A and 
B atoms in the alloy will not be spatially uniform. Some 
regions may contain significantly more A atoms than the 

1 Strictly, one needs to distinguish the bare distance from criti- 
cality that appears in J2j from the renormalized one that measures 
the distance from the true critical point. We suppress this difference 
because it is unimportant for our purposes. 



average. If the local A-concentration is sufficiently high, 
such regions will be locally in the magnetic phase even 
if the bulk system is nonmagnetic. Because the mag- 
netic fluctuations are overdamped, the quantum dynam- 
ics of sufficiently large such locally magnetic spatial re- 
gions completely freezes (for Ising symmetry |12j). At zero 
temperature, these rare regions thus develop static mag- 
netic order independently of each other. This destroys 
the sharp QPT by smearing [H1H31H3] and is manifest in 
a pronounced tail in the zero-temperature magnetization- 
composition curve [15] . 

At any nonzero temperature, the static magnetic order 
on individual, independent rare regions is destroyed be- 
cause they can fluctuate via thermal excitations. There- 
fore, a finite interaction between the rare regions of the 
order of the thermal energy is necessary to align them. 
This restores a conventional sharp phase transition at any 
nonzero temperature. However, the smeared character of 
the underlying QPT leads an unusual concentration de- 
pendence of the critical temperature T c which displays a 
tail towards large x [T31[T5] . 

The effects of disorder correlations can be easily under- 
stood at a qualitative level. For positive correlations, like 
atoms tend to cluster. This increases, at fixed composi- 
tion, the probability of finding large A-rich regions com- 
pared to the uncorrelated case. The tail of magnetization- 
composition curve therefore becomes larger (see Fig. Q}. 
In contrast, like atoms repel each other in the case of neg- 
ative correlations (anti-correlations). This decreases the 
probability of finding large A-rich regions and thus sup- 
presses the tail. 

Optimal fluctuation theory. To quantify the in- 
fluence of the disorder correlations, we now develop an 
optimal fluctuation theory [LUfTTj] . We focus on the "tail" 
of the smeared transition (large x) where a few rare re- 
gions show magnetic order but their interactions are weak 
because they are far apart. 

We roughly estimate the transition point in the alloy 
Ai—xB x , by setting the average distance from criticality to 
zero, r av = (I — x)r^ + xr& = 0. This defines the critical 
composition in "average-potential" approximation, 

x° c = -TA/(rB - va) ■ (3) 

For compositions x > x° c . static magnetic order can only 
develop on rare, atypical spatial regions with a higher than 
average A-concentration. Specifically, a single A-rich rare 
region of linear size Lrr can show magnetic order, if the 
local concentration x\ oc of B atoms is below some critical 
value x c . Because the rare region has a finite size, the 
critical concentration is shifted from the bulk value 
According to finite-size scaling [T61IT7] 

x c (L RR ) = x° c -DL^, (4) 

where 4> is the finite-size shift exponent and D is a 
non- universal constant. In a three-dimensional itinerant 
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magnet, <p takes the mean- field value of 2 because the 
clean transition is above its upper critical dimension. As 
x c {Lb.r) must be positive, a rare region must be larger 
than £ m i n = (D/x®) 1 ^ to show magnetic order. 

In the tail of the smeared transition, the magnetically 
ordered rare regions are far apart and interact only weakly 
To find the total magnetization M one can thus simply 
sum over all magnetically ordered rare regions. This gives 



M 



dL RR I dx\ oc P(N,X\ oc )m{N,X\ oc 



(5) 

where P(N,xi oc ) is the probability for finding a region of 
N sites and local composition xi oc (i-e--, a region contain- 
ing Nb = Nxioc atoms of type B), and m(N,xi oc ) is its 
magnetization 

Let us analyze the spatial distribution of atoms in the 
sample to determine the probability P(N, x\ oc ). Specifi- 
cally, let us assume that the random positions of the A and 
B atoms are positively correlated such that like atoms form 
clusters of typical correlation volume (number of lattice 
sites) Vdis «1 + a£^ is where £di S is the disorder correlation 
length and a is a geometric prefactor. The probabilities 
for finding A and B clusters in the sample are 1 — x and x, 
respectively. The number n c \ of correlation clusters con- 
tained in a large spatial region of N sites (N 3> Vdis) is 
approximately 



n c i « N/V dis = N/(l + a^ is ) 



(6) 



The probability P(N, xi oc ) for finding a region of N sites 
and local composition x\ oc is therefore equal to the prob- 
ability P c ius(^ci, »ib) for finding tib = xn c \ clusters of B 
atoms among all the n c \ clusters contained in the region. 
It can be modeled by a binomial distribution 



P c \u S {n Q \,n B ) = 

\n B J 



(7) 



We now distinguish two cases, (i) the regime where x is 
not much larger than x\ , and (ii) the far tail of transition 
at x — > 1. 

(i) If x is just slightly larger than rare regions are 
large and the probability ([7]) can be approximated by a 
Gaussian 



P 



1 



clus 



x /2nx(l - x)/n cl 



exp 



-n c v 



(•^loc 2/) 

2x{\ - x) 



(8) 



We estimate the integral (|S|) in saddle point approx- 
imation. Neglecting subleading contributions from 
m(N, xioc), we find that rare regions of size L^ R = D(2cj)~- 
d)/[d(x — x°)} and composition x c (Lrr) dominate the in- 
tegral. The resulting M(x) dependence reads 



M 



exp 



C (x- x° c ) 



0\2- 



(1 + < is ) *(1 - x) 



(9) 



where C = 2{D / d) d l ( >'{2(j) - d) d ^- 2 4> 2 is a non-universal 
constant. In this regime, varying the disorder correlation 



length thus modifies the non-universal prefactor of the ex- 
ponential dependence of M on x. 

(ii) An even more striking effect occurs in the tail of the 
transition for x — > 1. As rare regions cannot be large in this 
regime, the binomial distribution cannot be approxi- 
mated by a Gaussian. However, within saddle point ap- 
proximation, the integral ([5]) is dominated by rare regions 
containing only A atoms and having the minimum size per- 
mitting local order. Inserting L RR = L m ; n = (D/x®) 1 ^ 
and xioc = into ([S]), we find that the composition depen- 
dence of the magnetization is given by the power law, 



M ~ (1 - xf (x ->• 1) , 



(10) 



with f3 = aL^ lin /(l + a£dis)- In this regime, the disorder 
correlations thus modify the seeming critical exponent of 
the order parameter. The exponent value is given by the 
minimum number of correlation clusters necessary to form 
a magnetically ordered rare region. The results for uncor- 
related disorder |15) are recovered by substituting £dis = 
into © and (fTU|). 

So far we have assumed that a typical disorder correla- 
tion cluster of A atoms is smaller than the minimum rare 
region size required for magnetic order. For larger disorder 
correlation length £dis > L m i n , a single correlation cluster 
is already large enough to order magnetically. As a result, 
(almost) all A atoms contribute to the total magnetiza- 
tion. Correspondingly, the composition dependence of the 
order parameter is given by 



M ~ (1 - x) . 



(11) 



To combine the power laws (flO]) and (fTTj) for different 
ranges of £dis, we construct the heuristic formula 



f3 = (aLi in + < is )/(l + <D 



(12) 



which can be used to fit experimental data or simulation 
results. 

Other observables such as the finite-temperature phase 
boundary can be found in similar fashion. As discussed 
above, at T =/= 0, individual rare regions do not develop 
a static magnetization. Instead, global magnetic order 
arises via a conventional (sharp) phase transition at some 
transition temperature T c which can be estimated from 
the condition that the interaction energy between the rare 
regions is of the order of the thermal energy. To deter- 
mine the interaction energy, we note that in a metallic 
magnet, the rare-regions are coupled by an RKKY inter- 
action which falls off as r~ d with distance r. As the 
typical distance between neighboring rare regions behaves 
as r ~ M~ x / d [13], the composition dependence of the 
critical temperature is analogous to that of the magneti- 
zation. In particular, 

T c {x)~{l- X y (13) 

in the tail of the smeared transition, x — > 1. 
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Simulations. — We now verify and illustrate the the- 
oretical predictions by performing computer simulations 
of a toy model [I3].ll9j. Its Hamiltonian is motivated by 
the so-called quantum-to-classical mapping |18j which re- 
lates a quantum phase transitions in d space dimensions 
to a classical transition in d + 1 dimensions. The ex- 
tra space dimension corresponds to imaginary time in the 
quantum problem. Consequently, we consider a (3+1)- 
dimcnsional classical Ising model on a hypercubic lattice 
with three space dimensions and a single imaginary time- 
like dimension. The interaction in the time-like direction 
is long-ranged as the \oj n \ frequency dependence in @ cor- 
responds to a 1/r 2 in imaginary time. In the toy model, 
we replace this interaction by an infinite-range interaction 
in time direction, both on the same site and between spa- 
tial neighbors!! This correctly reproduces the smeared 
character of the phase transition due to static magnetic 
order on the rare regions. The Hamiltonian of the toy 
model takes the form 

H = — — JoSy iT S z , T i — — JySy tT Sy >T > , 

T (y,z),T,T' T y,T,T' 

(14) 

where y and z are space coordinates, r is the time-like 
coordinate, and S y , T = ±1. L T is the system size in time 
and (y, z) denotes pairs of nearest neighbors in space. J y 
is a binary random variable whose value, Jh or J;, is deter- 
mined by the type of atom on lattice site y. The values at 
different sites y and z are not independent, they are corre- 
lated according to some correlation function C(y — z) . The 
average concentrations of J^-sites and Jj-sites are 1—x and 
x, respectively. 

Treating the time-like dimension within mean-field the- 
ory, which is exact because of the infinite range of the 
interactions, a set of coupled nonlinear equations emerge 
for the local magnetizations m y = (l/L T ) Sy iT , 

m y = tanh — (J y m y + J2 z Jom z + h) . (15) 

J-cl 

Here, the z-sum is over the nearest neighbors of site y, 
and h is a tiny symmetry-breaking magnetic field. Ac- 
cording to the quantum-to-classical mapping, the classical 
temperature T c i is not related to the physical tempera- 
ture of the underlying quantum system (which is encoded 
in L T ) but rather some quantum control parameter that 
tunes the distance from the quantum phase transition. 

The local mean-field equations (fTBf can be solved effi- 
ciently in a self-consistency cycle. In the two clean limits 
with either J y = Jh or J y = J; for all y, the phase transi- 
tion occurs at T), = J/j+6 Jo and 2] = J;+6 Jo, respectively. 
We choose a classical temperature between Th and T\ and 
control the transition by changing the composition x. 

2 Even though the bare action ([TJ [2J does not have an interaction 
between spatial neighbors at different imaginary times r, such a cou- 
pling will be generated in perturbation theory (or under RG) from 
the short-range spatial interaction and the long-range interaction in 
time. 




Fig. 2: (Color online) Examples of the atom distribution in a 
plane of 256 2 sites for several values of the disorder correlation 
length £ dis = 0, 1.0, 2.0 from left to right (x = 0.5). 

To generate the correlated binary random variables rep- 
resenting the site occupations, a version of the Fourier- 
filtering method [20] is implemented. This method starts 
from uncorrclatcd Gaussian random numbers u y and turns 
them into correlated Gaussian random numbers v y charac- 
terized by some correlation function C(r). This is achieved 
by transforming the Fourier components u q of the uncor- 
rected random numbers according to 

5 q = [C(q)]*fiq, (16) 

where C(q) is the Fourier transform of C(r). The v y then 
undergo binary projection to determine the occupation of 
site y; the site is occupied by atom A if v y is greater than 
a composition-dependent threshold and by atom B if v y is 
less than the threshold. 

In the majority of our calculations, we focus on attrac- 
tive short-range disorder correlations of the form C(r) = 
exp (— f 2 /2£| is ). Figure [2] shows examples of the resulting 
atom distributions for several values of the disorder corre- 
lation length £dis- The formation of clusters of like atoms 
is clearly visible. 

We now discuss the results of the mean-field equations 
(]15[). Figure[3]presents the total magnetization M as func- 
tion of composition x for several values of £dis with all 
other parameters held constant. At a given composition x, 
the magnetization M increases significantly even for small 
£dis of the order of the lattice constant. Moreover, the 
seeming transition point (at which M appears to reach 0) 
rapidly moves towards larger compositions, almost reach- 
ing x = 1 for a correlation length £dis = 2. Inset (a) of 
Fig. [3] shows a plot of logM versus log(l — x) confirming 
the power-law behavior (fTOj) in the tail of the transition. 
The dependence on £dis of the exponents /3 extracted from 
these power laws is analyzed in inset (b) of Fig. [3l It can 
be fitted well with the heuristic formula. (|T2"|). 

In addition to the attractive (positive) correlations, we 
now briefly consider the case of anti-correlations (like 
atoms repel each other). We model the anti-correlations 
by a correlation function having values C(0) = 1, C(r) = 
— c for nearest neighbors, and C(r) = otherwise. The 
positive constant c controls the strength of the anti- 
correlations. A characteristic magnetization-composition 
curve for such anti-correlated disorder (with c = 1/6) is 
included in Fig. [3] The data show that the magnetization 
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Concentration, x 

Fig. 3: (Color online) Magnetization M vs. composition x for 
several values of the disorder correlation length £dis using one 
disorder realization of 256 3 sites, Ju = 20, J\ = 8, Jo = 1, 
T c i = 24.25, and h = 1CP 10 . Also shown is one curve for 
the case of anti-correlations (128 3 sites), for details see text. 
Inset (a): log-log plot of M vs. (1 — a;) confirming the power- 
law behavior in the tail of the smeared transition. The tail 
exponent /3 shown in inset (b) agrees very well with (|12[) as 
shown by the solid fit line. 

is reduced compared to the uncorrelated case, and the tail 
becomes less pronounced. Analogous simulations using 
different values of c show that this effect increases with 
increasing strength of the anti-correlations, as indicated 
in Fig. CD 

Conclusions. — In summary, we have studied the 
effects of spatially correlated disorder on smeared phase 
transitions. We have found that even short-range disor- 
der correlations (extending over just a few lattice con- 
stants) lead to qualitative modifications of the behavior 
at smeared transitions compared to the uncorrelated case, 
including changes in the exponents that characterize the 
order parameter and the critical temperature. In other 
words, systems with uncorrelated disorder and with short- 
range correlated disorder behave differently 

This is in marked contrast to critical points, at which 
uncorrelated disorder and short-range correlated disorder 
lead to the same critical behavior. (Long-range correla- 
tions do change the critical behavior [51[2"T].) What causes 
this difference between critical points and smeared tran- 
sitions? The reason is that critical behavior emerges in 
the limit of infinitely large length scales while smeared 
transitions are governed by a finite length scale, viz., the 
minimum size of ordered rare regions. This renders the 
renormalization group arguments underlying the general- 
ized Harris criterion [7J[8] inapplicable. 

The majority of our calculations are for the case of like 
atoms attracting each other. For these positive correla- 




Fig. 4: (Color online) Experimental temperature-composition 
phase diagrams of Sri-ajCa^RuO;!. Data from Hosaka et al. 
[23| . Schneider et al. [24], Wissinger et al. [25], and Khalifah et 
al. [26] are for thin films while those of Kiyama et al. [27], and 
Cao et al. [3H] are for bulk samples. Published magnetization 
curves show similar variations. 



tions, large locally ordered rare regions can form more 
easily than in the uncorrelated case. Thus, the tail of the 
smeared transition is enhanced; and the phase boundary 
as well as the magnetization curve move toward larger x 
as indicated in Fig. [T] We have also briefly considered 
the case of like atoms repulsing each other. These anti- 
correlations suppress the formation of large locally ordered 
rare regions compared to the uncorrelated case. As a re- 
sult, the phase boundary and the magnetization curve will 
move toward smaller x. In addition to short-range corre- 
lations, we have also studied long-range power-law corre- 
lations which arc interesting because they lead to a broad 
spectrum of cluster sizes. Detailed results will be pub- 
lished elsewhere [22] . 

Turning to experiment, our results imply that smeared 
phase transitions are very sensitive to slight short-range 
correlations in the spatial positions of impurities or de- 
fects. In particular, an analysis of the data in terms 
of critical exponents will give values that depend on 
these correlations. We believe that a possible realiza- 
tion of the effects discussed in this paper can be found 
in Sri_ K Caa;Ru03. This well-studied material undergoes 
a ferromagnetic QPT as a function of Ca concentration. 
Because Sri-xCa^RuOs is a metallic system with Ising 
spin symmetry, the transition is expected to be smeared 
[13] . Interestingly, the reported experimental phase dia- 
grams (see Fig. [4]) and magnetization curves show unusu- 
ally large variations. Not only does the apparent critical 
composition change between x ~ 0.5 and 1; the functional 
form of the magnetization curves also varies. Although 
part of these discrepancies may be due to the difference 
between film and bulk samples [25] , large variations within 
each sample type remain. We propose that disorder corre- 
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lations, i.e., clustering or anti-clustering of like atoms may 
be responsible for at least part of these variations. 

Finally, we emphasize that even though we have con- 
sidered the QPT in itinerant magnets as an example, our 
theory is very general and should be applicable to all phase 
transitions smeared by disorder including QPTs P9"H3"T] . 
classical transitions in layered systems (HO GUI an d non- 
equilibrium transitions |33j 

We thank I. Kezsmarki for helpful discussions. This 
work has been supported in part by the NSF under grant 
No. DMR- 0906566. 
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